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A  NEW  ADAPTIVE  ANTENNA  SYSTEM  FOR  COHERENT  SIGNALS  AND  INTERFERENCE 


T.  J.  Shan  and  T.  Kailatb 


Information  Systems  Lab..  Stanford  University 
Stanford.  CA  94305 


ABSTRACT 

In  this  paper  we  introduce  a  new  adap¬ 
tive  antenna  system  able  to  work  well 
even  when  the  desired  signal  and  the 
interference  are  coherent.  The  present 
adaptive  beamformers  fail  to  operate 
in  these  cases.  The  results  of  simula¬ 
tions  appear  to  confirm  the  theoretical 
predictions. 


I.  Introduction 

Since  the  pioneering  work  of  Howells[l]  . 
Applebaum[2]  and  Widrowl3]  .  there  has  been  con¬ 
siderable  activity  in  the  development  of  adaptive 
antenna  arrays  for  radar,  sonar,  communication, 
spectral  estimation,  etc. 

Though  the  details  differ  in  the  different  appli¬ 
cations.  the  main  assumptions  and  processing  algo¬ 
rithms  are  essentially  the  same.  In  particular,  a 
key  assumption  in  all  the  previously  cited  work  is 
that  the  interfering  signals  are  not  coherent  with 
(i.e..  do  not  have  fixed  phase  differences  from)  the 
desired  signal.  More  generally,  two  signal  will  be 
said  to  be  coherent  if  one  is  a  scaled  and  delayed 
replica  of  the  other.  Coherent  interference  can 
arise  when  multipath  propagation  is  present,  or 
when  "smart"  jammers  deliberately  induce 
coherent  interference,  e.g.  by  retrodirectlng  the 
signal  energy  to  the  receiver. 

Coherence  can  completely  destroy  the  perfor¬ 
mance  of  adaptive  array  systems.  We  shall  show 
this  by  computer  simulations  in  Section  IV. 

In  the  reference^]  Evans,  Johnson  and  Sun 
show  that  the  subaperture  sampling  or  spatial 
smoothing  idea,  as  they  call  it,  can  be  applied  to 
the  off-line  eigenstructure  based  method  (of  Bien- 
venu  and  of  Schmidt)  for  direction  finding.  This  is 
an  important  contribution,  the  main  idea  of  which 
we  independently  rediscovered  later  (see  Shan,  Wax 
and  Kailath[5]  ).  However,  adaptive  versions  of  this 
spatial  smoothing  scheme  were  not  obvious,  In  this 
paper  we  will  introduce  an  on-line  spatial  smooth¬ 
ing  algorithm  for  the  problem  of  adaptive  antenna 

TThil  work  woe  tup  ported  in  port  by  the  Air  Foree  Office  of  Scientific 
Reeeorch.  lir  Force  Syetomi  Commend  under  Contract  AF«<V6Z0-79-C- 
OOaa.  the  U.S.  Army  Reeeorch  Office,  under  Contract  DAAG29-7frC- 
0215,  and  by  tho  Joint  Service*  Program  ol  Stanford  University  under 
Contract  QAAC29-S1-K-00S7 
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array  system. 


II. 


The  Optimum  Weight  Vector  and  Beam 
Patterns  of  Beamforming  Array  for 
Noncoherent  and  Coherent  Incident  Signals 

It  is  easy  to  show  that  the  optimum  weight  vec¬ 
tor  lor  different  criteria  is  of  the  Iorm[6] 

«.p<  —  a(7o) 

where  Rn  =  Ex(t )x  (f)  is  the  covariance  matrix 
of  the  array  measurement  vector  *  ,  a  is  a  scalar 
constant,  and  n(T0)  is  so-called  'look''  direction 
vector. 

The  beam  pattern  of  a  array  of  M  sensors  is 
obtained  by  plotting 


where 


w.*„a(T) 


o(t)  =  (l  e-J"T  •  ■  t-ilU-Dvr-T 
t  =  d  cos  d/  c  ,  0  «  d  ts  2  w  . 


Explicit  analytical  expressions  are  difficult  to 
obtain,  but  we  can  obtain  insight  into  the  interfer¬ 
ence  rejection  properties  of  the  array  by  the  fol¬ 
lowing  asymptotic  (high  SNR)  analysis. 

We  shall  assume  that  we  are  interested  in  the 
signal  s(f)  in  the  known  look  direction  dg,  and 
that  this  signal  is  statistically  independent  of  the 
(interfering)  signals  from  the  other  K-l  unknown 
directions  . **-if 

U  in  addition  there  are  no  fixed  relations 
between  the  phases  of  signals,  or  more  generally  if 
none  of  the  signals  is  a  scaled  and  shifted  version  of 
any  other  signals  we  shall  say  that  the  signals  are 
completely  noncoherent  . 

Noncoherent  Signals 

Under  the  assumptions  of  statistical  indepen¬ 
dence  and  of  noncohercnce,  we  shall  express  the 
■irray  measurement  data  vector  as 

x(f)  =  a(r0)s(f)  ♦  Kj(f)  +  »(t) 
where  A  specifies  the  interference  directions 
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*  =  [»(T.)  •  •  •  •(T*-,)] 

and  s(0  .  the  array  input  signal  vector,  is  consist¬ 
ing  o!  desired  signal  s(t )  and  jamming  signals  j  (f ) 


*jri«** 

r  i*K 


Pi  =  e,»(T0) 


fs(Ol 

f 

■(0  =  = 

u(Ol 

The  covariance  of  x(<)  can  now  be  written 
R»  -  Po«(to)“’(to)  +  Rnn 

where 

Rnn  =  +  0*1 

and  we  have  assumed,  for  simplicity,  the  noise 
intensity,  cr*.  to  be  same  at  each  sensor. 

By  using  the  matrix  inversion  lemma,  we  can 
write 

wopi  =  <*R«a(r0) 

a  a|Rnn«(r0)  ~P§- ‘  ,(T°)R1-'a(r°)  1 
l  1  +  J>0*  (To)Rnn«(T0)  J 

*  P  Rnn*(To) 

where  p  is  a  scalar  constant. 

We  also  introduce  the  modal  representation 

*R i#**=  S\e4e;  . 


where  {X*  |  and  }e<  {  are  the  nonzero  eigenvalues 
and  the  corresponding  eigenvectors  of  the  M  x  M 
matrix  *Ra*',  which  will  have  rank  K-  t 
because  R }j  is  the  covariance  matrix  of  the  K  -  1 
noncoherent  signals. 

Finally  we  shall  also  assume  that  the  back¬ 
ground  measurement  noise  intensity  is  small  com¬ 
pared  to  the  signals  7(f).  so  that  we  shall  have 

_  \  »  o 2 


~»  — — 
oz  \  +  o* 


Then  we  can  write  * 
^ tpi  ~~  PRrtn  »(T0) 

iK-i  , 

-  a  I  I  


(K-\  ,  .  M  ,  J 

•  •I  K  *0  i-K  o* 


where  =  denotes  asymptotic  (  as  &*■* 0  )  equality. 
Now  by  construction,  the  direction  ve  .  -rs 
ia(T|),....a(Tjf_j)J  of  the  interfering  signals,  which 
are  the  columns  of  the  matrix  K.  lie  in  the  span  of 
the  first  K  -  1  eigenvectors  fe,.  .  .  .  and 

are  therefore  orthogonal  to  the  remaining  eigen¬ 
vectors  jejf . ejgj. 

Therefore,  we  shall  have 

w^a(r,)  il-f  Pie<*a(T,}  =  0  l  =  1 . K- 1  . 

Or'  imK 

so  that  the  beam  pattern  will  have  "deep  nulls"  in 
the  interference  directions.  In  the  look  direction, 
on  the  other  hand,  the  constraint  will  ensure  that 
we  have 

wop<«(T0)  =  1 

This  is  the  well  known  behavior  of  the  linear 
array,  which  can  be  approximately  achieved  by  a 
variety  of  adaptive  algorithms[6]  . 

However,  the  situation  deteriorates  badly  in 
the  coherent  case. 

Coherent  Signals 

If  the  signals  have  fixed  phase  differences, 
which  rally  means  equal  frequencies  and  fixed  p,. 
then  w-e  shall  have  the  representation 

As  =  a (r0)s(O  +  £  •(ri)ji(t) 


~  l«(To)  +  7i*(Tj)  +  •  •  +  ■)'Ar-j*(TA-_,)]s(t) 

where  the  (7tJ  are  fixed  complex  constants  given 
by 

7i  =  (Pi/Po)  e>(',,■r',  .  i  =  1 . K- 1 

In  this  case,  the  covariance  matrix,  AFsg'A’.  will 
have  rank  1.  so  that  it  will  have  only  one  nonzero 
eigenvalue  k,,  and  the  covariance  matrix  R„  will 
have  M  -  l  eigenvalues  equal  to  a*. 

Therefore  we  shall  have 
w.pi  =  aRx~'a(To) 

f  .  H  .  J 

=<*  T — 7 — Te>e»*  +  2  3*c>e»|  *Ct0) 

jk,  +  a4  «-z  r  ] 

=  *  £  -Vfe.afT,,))*!, 

«-z  o 

All  we  can  say  here  is  that  the  linear  combination 


b  :=  a(r0)  +  +  •  •  •  + 


z<»(t)  =  [*,  ■  ■  wof 


will  be  orthogonal  to  the  (e2 . e#j.  This  does 

not,  however,  imply  that  the  same  will  be  true  of 

the  Ja(r,) . a(T*_,){  individually,  and  therefore 

there  will  not  in  general  be  any  nulls  in  the  direc¬ 
tions  of  the  interfering  signals. 

Due  to  the  constraint,  we  shall  have 
w.j>«»(To)  *  1  . 

but  this  is  of  small  comfort,  because  the  actual 
array  output  will  be 

V(0  =  w.V(X(t) 


Define 

R<*>  =  £Tz<*>z<*>' 

and  the  spatial  smoothed  correlation  matrix: 

h  =  ^-tRg’ 

Pkml 

Then  we  can  prove  (  see  Reference[7, 5]  )  that  R 
will  have  the  form 

R  =  ASA  '  +  a*  I 

where 

S  will  have  rank  K  if  and  only  if  p  a  K  . 


u 

where  we  recall  that  b  :=  a(r0)  +  £  7(a(T()  lies 

along  Cj  and  is  orthogonal  to  fe2 . 

Therefore  there  will  be  no  signal  output  from  the 
conventional  array  when  the  signals  are  coherent. 

A  Way  Out 

This  analysis  also  makes  clear  what  is  neces¬ 
sary  to  rescue  the  situation:  we  must  somehow 
restore  the  rank  of  the  covariance  matrix 
£(As)(As)’  to  being  K.  Then  the  noise-alone 
eigenvectors  will  be  orthogonal  to  all  the  vectors  in 
the  space  of  the  signals  (desired  signal  and 
interfering  signals)  and  the  beam  pattern  will  have 
nulls  in  the  directions  of  the  interfering  signals.  A 
simple  scheme  for  achieving  this  rank  restoration 
with  an  adaptive  algorithm  is  proposed  in  the  next 
section. 


HI.  A  New  Adaptive  Antenna  Array  System 

We  shall  describe  a  preprocessing  scheme  for 
the  sensor  outputs  that  will  restore  the  rank  of  the 
signal  covariance  matrix  to  K  even  if  the  signals 
are  completely  coherent  with  each  other. 

The  scheme  is  based  on  combining  measure¬ 
ments  from  overlapping  subarrays. 

Given  the  M  sensor  outputs  at  any  time 
instant, 

*(0  =  l*i(0  **(f)jr  . 

define  p  subsets  (recall  that  K  is  the  number  of 
sources) 

*(,,(0  =  l*i(0  •  *XM<‘)]r 


Once  S  has  rank  K,  then  the  noise  eigenvec¬ 
tors  will  be  orthogonal  to  the  columns  of  A  and  by 
the  analysis  of  Section  III,  will  give  nulls  in  the 
interference  directions.  The  definition  of  a1**  shows 
that  K  +  p  =  M  and  this  combines  with  the  con¬ 
straint  on  p  to  require  that 

M  =  K  +  p  > 

Therefore  for  this  scheme  to  work,  we  must  have  at 
least  twice  as  many  sensors  as  signal  sources  in 
this  case. 

Coherent  Subgroups 

If  there  are  some  coherent  source  inputs  and 
some  noncoherent  inputs,  we  should  divide  the 
sources  into  C  noncoherent  subgroups  within 
each  of  which  the  inputs  are  completely  coherent. 
Then  we  form  a  matrix 

R  —  Rj  +  +  Rc 


To  destroy  the  coherency  in  all  of  the  groups,  the 
total  number  of  subgroups  must  be  at  least  equal 
to  the  size  of  the  largest  subgroup  (see  refer- 
encel?]  ). 

With  these  results  in  hand,  we  can  now  explain 
how  to  do  adaptive  processing  of  the  sensor  out¬ 
puts. 


Adaptive  Processing 

'Ve  can  rewrite  the  expression  for  the  estimate 
of  If  from  N  data  snapshots  as 

f 


=  -£  *;*>*;*>• 


*(*>(0  =  l*z(0  WOjr 


-hi  i 

J«1*M  1  3 


ihere  the  subscript  }  denotes  the  j-th.  lime 


instant,  i.e. 


=  [**(/) . **  0)] 


and  the  ’hat’  signifies  an  estimated  quantity. 
Then 

Ratm  =  t  *J*,*/*>' 

(iv  +  UP  jmlkml 


-  mhwU^"' 

=  l.w  i; Rv  *  ,V.1>S’M„.W, 


This  expression  suggests  that  we  can  recursively 
update  the  inverse  of  R  by  using  the  matrix  inver¬ 
sion  lemma  iteratively  p  times  once  for  each 
It  follows  that  we  can  also  use  approximate 
gradient-type  adaptive  algorithms  (e.g.  the  LMS 
algorithm  of  Widrow  and  HofT)  to  update  the 
weights  of  the  adaptive  processor. 

Fig.  la  shows  how  to  form  the  spatial  data  sub¬ 
set  from  one  ’snapshot’.  A  flow  diagram  of  the  pro¬ 
cedure  is  shown  in  Fig.  lb.  At  each  time  instant,  the 
’snapshot’  of  M  data  samples  is  divided  into  over¬ 
lapping  subgroups  of  K  samples  each:  these  sub¬ 
groups  are  then  fed  in  succession  into  the  adaptive 
processor,  which  updates  a  K-dimensional  (if  each 
sensor  is  followed  by  a  N  tapped-delay-line.  then  a 
KxN  -dimensional)  weight  vector  each  time.  After 
all  the  subgroups  have  been  processed,  the  same 
procedure  is  repeated  with  the  next  data 
’snapshot’. 

It  is  important  to  note  that  the  processor  can 
be  any  kind  of  adaptive  processor  using  any  adap¬ 
tive  algorithm  and  any  array  structure,  e.g.  the 
Howells-Applebaum  or  the  Frost  arrays  could  be 
used. 


IV.  The  Computer  Simulation  Results 

To  examine  the  performance  of  the  suggested 
beamformers  in  coherent  receiving  environments, 
several  computer  simulations  have  been  carried 
out,  with  that  confirm  the  theoretical  predictions. 

In  our  first  example,  we  have  a  signal 
0.1sin0.2nt  arriving  at  90*  and  four  coherent 
interfering  signals,  I0sin0.2nf,  arriving  at  25'  , 
75', 120®  and  160®  .  Fig. 2a  shows  the  beampattern 
of  a  conventional  Howells-Applebaum  array  with 
six  sensors.  Our  new  scheme  uses  10  sensors  with 
subarrays  of  size  six,  adapted  with  the  same 
Widrow-Hoff  LMS  algorithm,  and  gives  the  beampat¬ 
tern  shown  in  Fig. 2b.  The  input  signal  waveform  (at 
arrival  angle  90*  )is  shown  in  Fig. 3a.  with  the  out¬ 
put  of  the  conventional  array  in  Fig. 3b  and  that  of 
our  new  system  in  Fig. 3c  -  notice  the  big  difference 
in  the  scales  of  these  figures. 


In  the  second  example,  we  consider  the  case  of 
a  wide-band  signal  of  unit  power  level  with  central 
frequency  0  25,  bandwidth  0.1,  and  a  narrow-band 
interfering  signal  with  frequency  0.25.  at  power 
level  10  arriving  at  45®.  A  conventional  Frost  array 
will  suffer  considerable  spectral  and  waveform  dis¬ 
tortion  (signal  cancellation)  in  this  case;  again  the 
new  adaptive  algorithm  gives  much  better  perfor¬ 
mance.  The  input  signal  spectrum  is  shown  in 
Fig  . 4a;  Fig. 4b  and  Fig.4c  show  the  output  spectra  of 
the  Frost  array  and  of  the  new  scheme. 

Several  other  simulations  have  shown  similar 
results. 


V.  Conclusions 

Conventional  adaptive  antenna  arrays  perform 
very  poorly  in  coherent  receiving  environments.  If 
the  received  signal  is  coherent  with  one  interfer¬ 
ence,  the  signal  will  be  canceled  out  on  the  output 
of  the  adaptive  antenna  system,  which  will  there¬ 
fore  totally  fail  to  operate  as  a  receiving  unit.  The 
suggested  new  adaptive  beamforming  system  is 
able  to  overcome  this  degradation  of  performance 
in  coherent  receiving  environments,  without  con¬ 
siderably  increasing  the  complexity  of  the  system 
structure  or  the  computational  burden.  The  new 
array  structure  can  be  applied  in  conjunction  with 
any  of  the  adaptive  algorithms  and  structures  of 
current  adaptive  arrays,  and  successfully 
separates  the  coherent  array  inputs,  as  shown  by 
theoretical  analysis  and  simulation  results. 
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Figure  la:  Subaperture  sampling. 


Figure  lb:  Flow  diagram  of  the  new  scheme. 
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Figure  3b:  Outputs  signals  of  the 
Howells-Applebaum  Array. 
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Figure  3c:  Output  signal  of  the  new  array. 
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Figure  2a $  Beampattern  of  the  Hovel Is-Applebau3 
array  with  coherent  Inputs. 


Figure  2b:  Beampattern  of  the  new  array  with 
the  sam  v*  inputs. 
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Figure  4a:  Input  signal  spectrum. 
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Figure  4b:  Frost  array  output  spectrum. 
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Figure  4c:  New  array  output  spectrum. 
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